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ABSTRACT 

We have undertaken a search for millimetre-wave band absorption (through the CO and 
HCO + rotational transitions) in the host galaxies of reddened radio sources (z = 0.405 — 
1.802). Despite the colour selection (optical-near infrared colours of V — K > 5 in all but 
one source), no absorption was found in any of the eight quasars for which the background 
continuum flux was detected. On the basis of the previous (mostly intervening) H2 and OH 
detections, the limits reached here and in some previous surveys should be deep enough to 
' detect molecular absorption according to their V — K colours. However, our survey makes 

the assumption that the reddening is associated with dust close to the emission redshift of 
the quasar and that the narrow millimetre component of this emission is intercepted by the 
compact molecular cores. By using the known millimetre absorbers to define the colour depth 
and comparing this with the ultra-violet luminosities of the sources, we find that, even //"these 
assumptions are valid, only twelve of the forty objects (mainly from this work) are poten- 
tially detectable. This is assuming an excitation temperature of T x = 10 K at z = 0, with the 
number decreasing with increasing temperatures (to zero detectable at T x > 100 K). 
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1 INTRODUCTION 

Millimetre-wave observations of molecular absorption systems 
along the sight-lines to distant quasars provide a powerful probe 
of the cold, dense, star forming gas in the distant Universe. Fur- 
thermore, comparison of the redshifts of the rotational transitions 
of the molecules with those of the spin-flip transition of H I, as well 
as the electronic optical/UV transitions of metal ions, can be used 
to determine high redshift values of the fundamental constants, to 
at le ast an order of magn itude the sensitivity of purely optical data 
(see lCurran et al. I l2004al) . However, des pite much searching, only 
four su c h syste ms are currently known (IWiklind & Combesll 19951 
1 1996al lbl lT997h . the highest redshift being at 2 abs = 0.89. Of these, 
two are intervening systems (gravitational lenses towards more dis- 
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tant quasars), with the other two systems arising through absorption 
within the host galaxy of the quasar. Subsequent searches at the red- 
shifts of known high column density H I absorption systems, inter- 
vening the sight-lines to more distant quasi-stellar objects (QSOs), 
have also failed to detect molecular absoiption in the millimetre- 
band |Cumn|et^|20043 and references therein), despite the possi- 
bility that these so-called damped Lyman-a systems (DLAsjj may 
account f or more than 80 per-ce nt of the neutral gas content in the 
Universe jProchaska et al.ll2005l) . 

DLAs are, however, not devoid of molecular gas: To date, the 
Lyman and Werner ultra-violet bands of H2 have been detected in 



1 These have neutral hydrogen column densities of A r Hi > 2 X 10 20 cm -2 
and are usually detected at 2 abs ~ 1.8, where the Lyman-a transition is 
redshifted in to the optical band. 
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Figure 1. The H2 column density (top) and the molecular fraction (middle) 
versus the observed frame optical-near infrared colour (where available) for 
high redshift molecular absorption systems. The circles represent the H2- 
bearing DLAs (all optically selected intervening absorbers) and the squares 
and stars the OH absorbers (radio selected), w ith the inset in the middle 
panel showing the normalised OH line strength l lCurran et af1l2006l) . These 
are comprised of the four systems originally identified in millimetre-wave 
transitions (with the least-squares fit to these shown) plus the gravitational 
lens at z^ s = 0.764 towards 0132-097 (detected in OH decimetre but not 
HCO + millimetre absorption. IKanekar et al .120051) . P(t) shows Kendall's 
r two-sided probability of the observed distribution occuring by chance 
and S(t) the significance of this assuming Gaussi an statistics. The bot- 
tom panel shows the distribution of SPSS quasars ^Schneider et alj|2007l 
after conversion to V magnitudes, iFukugita et al.l Il996h where the filled 
histogram shows those which have been detected in the Very Large Array's 
"Faint Images of th e Radio Sky at Twenty Centimetres" (FIRST) survey 
lBeckeretafll995t) . 



decimetre band absorbers have molecular fractions T ~ 0.6 — 1 
and optical-near-infrared colours of V — K > 5. 

The correlations in Fig. [7] present strong evidence that the 
quasar light is reddened by dust in the foreground absorber: Since 
the presence of the dust is necessary to prevent the dissociation 
of the molecular gas by the ambient ultra-violet field, the molec- 
ular fraction is expected to be correlated with the dust abundance, 
as observed. The paucity of millimetre-wave band absorption can 
therefore be attributed to the traditional optical selection of targets 
biasing towards absorbers of low dust content and therefore low 
molecular fractions. 

The fact that intervening absorbers are usually found through 
optical spectroscopy, yielding a redshift but also giving the above 
bias against dusty objects, means that millimetre-band searches 
of known intervening a bsorbers have generally been unsuccess- 
ful dCurran et al 1l2004bl and references therein). An alternative tar- 
get for molecular absorption is towards the fainter "red quasars", 
where the red colour may indicate an intervening column of dust. 
However, due to the relatively narrow bandwidths in the millimetre 
band (see Sect. I4.2t . such an approac h is currently only p ractical 
at longer (decimetre) wavelengths (see lCurran et al.llioilbl) . In the 
absence of any known intervening absorbers, selecting the quasar 
itself gives a redshift (z cm ) to which to tune the receiver. Naturally, 
such a selection of targets prevents any useful comparison with the 
optical redshifts, in order to measure the values of the fundamental 
constants, although any detections could be followed up in 21 -cm, 
giving the redshift of the spin-flip transiti on of H I. 

In lCurran etaLI d2006L I2OO8L 1201 1 ah . we presented the results 
from our decimetre-wave searches for such "associated" (OH and 
Hi) absorption and here present the results of our millimetre- wave 
survey for associated absorption. 



2 OBSERVATIONS 
2.1 Target selection 

As per ICurran et al. I d2006l) . our sources were selected 
from the Parke s Hal f-Jansky Flat-spectrum Sample (PHFS, 
iDrinkwater et al.l [l99 7 f\ on th e basis of their optical-near-IR 
photometry dFrancis et al. I l2000l) . From these, we selected the 30 
reddest sources (which correspond to an extinction of Ay w 4.1), 
in which the emission redshift of the quasar (z em ) would place a 
strong absorption line (CO or HCO + ) into the 3-mm band. After 
culling these further, by selecting those of S > —30° (thus being 
observable from northern latitudes jf| and with 3-mm flux densities 
expected to be > 100 mJy, the ten objects listed in Table 13.11 
remained. 



19 DLAs (seelNoterdaeme et alfeOOfcTl Ijorgenson etal]|2009l and 



ISrianand et al.1 2010h . These, however, have molecular abundances 
which are generally much lower t han those detectable with current 
microwave and radio telescopes (Curran et al.l l2004bl and Fig. Q] 
top). Furthermore, in lCurran et al. I2OO6I) we showed that the H2- 
bearing DLAs have molecular fractions of T = tttt — t%t — ~ 

2« H2 +« HI 

10~ 7 — 0.3 and V — K < 4 (Fig.[T] middle), i.e. in the same range 
as a "typical" QSO (Fig.[Q bottomfl whereas the millimetre and 



2 On e of which, J1439+1 1 3, has also been detected in the CO A — X UV 
band ISrianand et al.l2008t) . 

3 V — K = 2.88 ± 1.04 in general and 3.05 ± 0.97 if radio-loud. 



2.2 The IRAM 30-m observations 

From December 2003 to February 2004 we observed three of 
the targets with the IRAM (Institut de Radio Astronomie Mil- 
limetrique) 30-m telescope at Pico Veleta, Spain. We used two 3- 
mm SIS receivers (A100 and B100), tuned to the redshifted fre- 
quencies of the molecular transitions (see Table 13. p . The obser- 



4 With the additio n of 0500+019, in cluded since it has been detected in 
21-cm absorption l lCarilli et al.lll998l) . We also inc luded J0906+4952 & 
J1341+3301, which are two very red sources from iGlikman et all i2004l) 
[Sect. [23]. 

5 WemissSEST. 
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vations were done with a nutating subreflector, switching symmet- 
rically ±90" in azimuth with a frequency of 0.5 Hz. The contin- 
uum levels of the observed sources were determined using a con- 
tinuum backend and increasing the subreflector switch frequency 
to 2 Hz. The image sideband rejection of the receivers were high, 
of the order of 20 dB (single sideband). System temperatures typ- 
ically ranged from 120 K (90 GHz) to 180 K (1 13 GHz). We used 
the full 1 GHz backend with broad (1MHz) filterbanks and nar- 
row band autocorrelators, the former giving a channel spacing of 
« 3 - 4 km s _1 (TableEm over a bandwidth of « 1500 krns" 1 , 
which should be sufficient to cover uncertainties in the emission 
redshifts, all of which are known to at least three decimal places 
(Az = ±0.001 corresponds to Av « ±100 - 200 km s _1 for 
our sample). The pointing of the telescope was checked regularly 
on nearby continuum sources. Typical pointing corrections were 
5 — 10". The focus was checked regularly on Mars and Saturn. The 
Half Power Beam Width (HPBW) at 95 GHz is 26" . 

We prioritised the three targets according to the 3-mm flux 
densities estimated from an interpolation of the decimetre and near- 
infrared values: 12.8 hours of integration on 0500+019 (5 es t ~ 
0.32 Jy), 13.6 hours on 1430-155 (S cst £ 0.08 Jy) and 7.2 hours 
on 1504-166 (S DS t ~ 0.8 Jy), which, upon comparison with the 
observed values (Table [3~Q , were reasonable estimates. The data 
were reduced with the GILDA^fl software package. 



2.3 The NRO 45-m observations 

The remainder of the sample was observed with the Nobeyama 
Radio Observatory's 45-m telescope in March 2004. We used the 
H28/32 (1-cm) and the S80, S100 (3-mm) receivers to observe 
the J = — ► 1 and 1 — > 2 transitions over a range of red- 
shifts (Table 13.11 ). The observations were performed in position 
switching mode, with an integration time of 20 s for each scan. 
The antenna temperature, , was obtained by the chopper-wheel 
method, which corrected for atmospheric and ohmic losses. System 
temperatures typically ranged from 200 K (1-cm) to 500 K (3-mm). 
We used the AOS backend over 250 MHz split over 2048 chan- 
nels, which gave channel spacings of 0.37 (3-mm) kms -1 to 1.1 
km s" 1 (1-cm) with bandwidths of w 700 and « 2000 km s" 1 , 
respectively. The pointing of the telescope was checked by observ- 
ing SiO tnaser sources with the H40 (40 GHz band) receiver and 
the corrections were ^ 10" . 

The 3-mm observations were performed when the weather 
was clearest, with total integration times of 2 hours for 0213- 
026, 5 hours for 0454+066, 3.7 hours for 1107-187 and 8.5 hours 
for 1706+006 (for which we could not determine the flux den- 
sity, Fig. During less than ideal weather conditions, we ob- 
served in the 1- cm band and in c luded the radio detected, optically 
dim sources of iGlikman et al.l J2004I) . where HCO+ ->• 1 is 
redshifted into this band — the FIRST-2MASS reddened quasars 
J0906+4952 (SDSSJ09065 1.49+495235.9) and FTM J 134 1+3301. 
The data were reduced with the NEWSTAR package and, like the 
IRAM results, upon the removal of a low order baseline and 
smoothing no absorption features were apparent in the spectra. 



3 RESULTS 



3.1 Observational results 



In Fig. [2] we show the reduced spectra and in Table |3~T| we sum- 
marise our observational results. The limit to the total column den- 
sity of each molecule is calculated from 
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where v is the rest frequency of the J — > J + 1 transition, gj+i and 
Aj+i^j are the statistical weight and the Einstein A-coefficient of 
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is the 



the transition, respectively, and Q = X/ /-o 9 J e 
partition functiorjj. In the case of the four kno wn systems, the cov- 
ering factor is expected t o be close to unity dWiklind & Combesl 
ll994Lll995lll996bL [ 19981) so, as in the optically thin regime, this 
can be written outside of the integral. However, unlike lower fre- 
quency searches (e.g. OH 18-cm, 2 n 3/ / 2 ) I in the millimetre regime 
even low excitation temperatures give fcT x ~ Ej ~ hu and so the 
column density cannot be approximated via a linear dependence 
on this. Therefore, unlike the decimetre se arches for associated ab- 
sorption dCurran et al . 2006, 2008, 20TTah, we must assume an ex- 
citation temperature. 

In Table [2] we apply Equation Q] to the velocity integrated op- 
tical depths given in the references in order to derive the column 
densities and excitation temperatures of the f our known systems. 
These are found to be close to those derived bv lWiklind & Combesl 
jl 99511 1996al lbl l 19971) and so for our targets we adopt an excitation 
temperature of T x = 10 K at z = 0. 

Although the excitation temperatures are similar, the four 
known millimetre absorbers exhibit a range of full-width half max- 
ima (FWHMs), ranging from FWHM ;s 2 kms -1 (1413+135, 
IWiklind & Combesll 19971) to ~ 70 km s 1 (t he main component 



towards 1504+377. IWiklind & Combeslll996bh . For these, plus the 
absorber towards 0132-097 (Sect. 1), it has been shown that the 
OH 18-cm ( 2 n, y 2 J = 3/2) FWH Ms are similar to those of the H I 
21-cm profiles dCurran et al .1120071) . although we cannot yet unam- 
biguously state that millimetre widths are correlated with those of 
the decimetre lines (Fig. [3}- In any case, of all of our targets, only 
on e (0500+019) has b een detected in 21-cm absoiption (see table 6 
of ICurran et al. I l2008l and Sect. 14.2) , and so we would have no in- 
formation on what the width of each of the undetected lines should 
be. If this were known, we would smooth the spectral resolution 
(Av) of the data to the FWHM of the line, thus giving the best 
possible optical depth limit, based upon the detection of absorption 
in a single channel. Multiplying the resulting optical depth, scales 
the velocity integrated optical depth to r x FWHM/ Av, giving an 
overall scaling of ^/FWHM/Au. 

Choosing either the minimum and maximum FWHM of the 
four known absorbers, therefore gives a possible range in scaling 
factors of a/ 70/2 « 6, and so the choice of either FWHM « 2 
kms -1 or w 70 kms -1 could lead to a significant over/under- 
estimate in the profile width. We therefore use the optical-near- 
infrared colour, which appears to be correlated with the profile 



6 http://www.iram.fr/IRAMFR/GILDAS/ 

7 The highest frequency flux measureme nt available for 1706+006 is 0.44 
Jy at 5.0 GHz IWright & Otrupceklll990l) . 



8 The Einstein A-coefficients are taken from IChandra et all dl99j 1 19961) 
or derived from the permanent electric dipo le moment of the m olecule, ob- 
tained from the JPL Spectral Line Catalog jPickett et all 19981) , along with 
the energy of each level, Ej . An on-line column density calculator based on 
Equationfjjis available from http://www.phys.unsw.edu.au/~sjc/column/ 
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Figure 2. The observed spectra before baseline removal. The ordinate in each spectrum shows the antenna temperature in Tj£ [K] and the abscissa the 
barycentric frequency [GHz], with the top axis showing the velocity offset from the central tuned frequency [kms -1 ]. All are shown at a resolution of 
4 km s" 1 , the coarsest of the observations (see Table l3.lV The negative spike in each of the IRAM spectra are due to bad channels (246-250) and the negative 
flux for 1706+006 is due to the long (20 s) integrations with the NRO giving errors in the continuum flux measurement. 



Table 2. The excitation temperatures and column densities derived from the listed references. In the case where the T x entiy is 
empty, the column density is calculated at the excitation temperature of the other molecule. 



System 


-^abs 


Reference 




CO 


HCO+ 








T x [K] 


iVtotal [cm~ 2 ] 


T x [K] AT total [cm- 2 ] 


0218+357 


0.68466 


Wiklind & Combes (1995) 


9 


6.3 + 0.3 x 10 16 


5.1 + 0.2X10 13 


1413+135 


0.24671 


Wiklind & Combes (1997) 




1.54 + 0.04 x 10 16 


8 2.4 + 0.4 x 10 13 


1504+377 


0.67335 


Wiklind & Combes (1996b) 


16 


5.9 ± 1.0 x 10 16 


13 6.5 + 0.1 x 10 13 


1830-211 


0.88582 


Wiklind & Combes (1996a, 1998) 




si 6 x 10 17 


8 2.9 + 0.3X10 14 
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Figure 3. The FWHM of the H I and OH profiles versus the optical-near- 
infrared colour for the five OH absorbers. The black symbols show the ap- 
proximate millimetre FWHM versus the optical-near-infrared colour, for 
the four of the five detected in millimetre lines,. The inset shows the mil- 
limetre detections only with the line showing the fit used to estimate the 
FWHM for our targets. 



width in the four known systems (Fig. [31 L in order to estimate 
what the FWHM of an undetected absorber should be. 

As seen from the figure, the associated absorber in 1504+377 
represents a stray point with a FWHM of ~ 70 km s _1 . We there- 
fore estimate the profile widths from the least-squares-fit to the 
other three systems (shown in the inset of Fig. |3}. The fact that 
this is based upon only three points, with the removal of one which 
we simply do not like, limits the robustness of this method, al- 
though it should give a better handle on the profile width than 
a simple assumption. We therefore use this fit, log 10 FWHM ~ 
0.86 (V-K)- 3.84 (for 2 sC FWHM sC 70 km s" 1 ), to estimate 
the column density limits (given in Table[3~llFI 



3.2 Comparisons with the known high redshift molecular 
absorbers 

In Fig. [4] we show the CO and HC O + column densities normalised 
by that of the Hi 21 -cm (from ICarilli et all Il992l Il993l Il998t 
IChengaluretaLlll999h against the optical-near-IR colour, for the 
four known millimetre absorbers. From this it is difficult to as- 
certain whether there is a correlation between the normalised 
millimetre-wave column density and the V — K colour, as i s found 
for the normalised OH column densities in ICurran et al. | ( 120061) 
[Fig-U). Furtherm ore, 21 -cm line stren gths are not available for the 
sample although in ICurranetal.ld2006h there was still a, albeit more 
scattered, correlation between the un-normalised OH 2 n 3( / 2 J r = 
3/2 column density and the V — K colour, which appears to apply 
to all of the molecular absorbers (Fig.Q] top). Therefore in Fig.|5]we 
show the un-normalised CO and HCO + column densities against 
the colour for our targets, as well as those from previous searches of 
associated molecular absorption dWiklind & Combeslll995lll996bl : 



9 This suggests that the path length through the dust responsible for the 
reddening may be correlated with the orientation of a rotating cloud com- 
plex, such as the disk of a galaxy (Curran & Whiting, in prep.). There is a 
1% probability of this correlation arising by chance (for H I, Fig. [5), i.e. a 
3.30ct significance assuming Gaussian statistics. 

10 In order to guard against over-compensating, we use the minimum and 
maximum observed profile widths of 2 and 70 km in cases where the 
estimated FWHMs are less or greater than these limiting values. 
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Figure 5. The CO (left) and HCO+ (right) column densities versus optical-near-IR colour. The symbols are as per Fig. [4] with the filled symbols show i ng the 
four known systems (the l i ne sho ws the least-squares fit to th ese) and the unfilled symbols showing the limits from the searches bv lWiklind & Combed Jl995l 
ll996bh;|Prinkwater et aTUl996l) fbluel iMurphv et alj J20031) (specifically 0727-115 at 2 cm = 1.591, see Sect. [4/2) [red] and the millimetre wave searches 
of ICurran et all 1 12006 ) [green] (actually HCN in 0335-122, since this is the only source for which V and K are available) and this paper [orange]. All of the 
limits are at the 3(7 level per channel recalculated for the expected FWHM, except when outside the range of the four known, where either 2 or 70 km s _1 is 
used (Sect. l3~TV Note that this range may add a potential uncertainty of 0.78 dex to the column density estimates (Sect. l3~71 although see footnotes ll2l and lT3l . 
in addition to further uncertainty from the assumption of T x = 10 K at z = (see Sect. 14. 3t . 
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Figure 4. The normalised (assuming T s /f = 100 K) CO (black) and 
HCO+ (coloured) column densities versus optical-near-IR colour. The dot- 
ted lines show the least-squares fit for each of the molecules with the dashed 
line showing the composite fit. 



Drinkwater et al ] fl996l: iMurohvetaLlEool ICurran et alEooOT 
From these, we see that many of the limits (especially the new ones 
reported in this paper, shown in orange), should be more than suffi- 
cient to detect molecular absorption, particularly HCO + . Incorpo- 
rating the upper limits o f the column dens ities, via the ASURV sur- 
vival analysis package dlsobe et alii 19861) , gives Kendall's r two- 
sided probabilities of P(r) = 0.0092 (n = 26) [i.e a signifi- 
cance of S(t) = 2.60a for COH and P(t) = 0.025 (n = 27) 



In the 



plots we give Iwiklind & Combed il995l Il996bl) ; 

iDrinkwater et all i ll 9961) the same colour as the latter paper includes 
some of the results of the first two. 



Bearing in mind that these limits are subject to the FWHM estimates 
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(Sect. I3.lt , we have also derived the probabilities when applying the two 
FWHM extrema of the detected absorbers, 2 and 70 km s —1 . Not generally 
being as red as the detected sample, the values do not change for FWHM= 
2 km s _1 , and for 70 km s _1 , P(r) 70 = 0.0075 & S(r) 70 = 2.68(7 
(n = 28). 



[S(t) = 2.24(7 for HCO+H of the column density-colour cor- 
relations occuri ng by chance. For the five known OH absorbers, 
S(r) = 1.96cr <Curranetalj2011al) . which falls to 1.36cr for those 
four detected in millimetre transitions^ So although the search re- 
sults do not lessen the significance of the column density-colour 
correlations (but actually increase them), many of the limits are be- 
low those expected for a detection and so we discuss other possible 
reasons as to why millimetre-wave spectral lines were undetected 
in this and the previous surveys. 



4 DISCUSSION: POSSIBLE EXPLANATIONS FOR 
THESE AND THE PREVIOUS NON-DETECTIONS 

4.1 Millimetre-wave covering factors 

While there is an apparent correlation between the normalised col- 
umn density and the optical-near infrared colour, Fig. Q] actually 
shows the ratio of the OH 18-cm velocity integrated optical depth 
to that of the Hi 21 -cm. Therefore in order for the fit to be accu- 
rate for the normalised column densities, the four known millimetre 
absorbers must also have similar spin and excitation temperatures, 
as well as H I and OH covering factors (or at least a similar ratio 
of these quantities). While the excitation temperatures can be well 
constrained from the millimetre transitions (Table |2] and the spin 
temperatures can be removed by not normalising the column densi- 
ties (Fig.[4]cf.[5}, the covering factor, /, is ingrained into the optical 
depth via, r = — ln(l — 5%). Since the CO and HCO + lines are 

optically thick in the four known systems, — In (l — -^J ^ y^, 
thus not allowing / to be taken outside of the integral in EquationQ] 
although, as noted in Sect. |3.1| if / ~ 1 this would have little effect. 

In any case, the correlation shown in Fig.|4]does not hold for 
all five of the absorbers: Towards the z e l n = 2.215 quasar PMN 
J0134-0931 (0132-097), Ikanekar et al] J20051) detected OH at a 



13 Again, applying FWHM= 2 km s 1 introduces no change to the cor- 
relation and P(t) 7u = 0.022 & 5(r) 7 o = 2.29ct (n = 30). 

14 S(t) = 2.25(7 for the five known OH absorbers plus the OH limits 
ICurran et ai1l201 id) . 
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column density of 7V OH = 5.6 x 10 14 . (T x /f) cm in the inter- 
vening z a bs = 0.765 gravitational lens, although HCO + 1—^2 
was undetected to r < 0.0 7. This gives a col umn density ra- 
tio of iVoH ~ 5600iV HCO+ . ICurran et ail d2007l) have suggested 
that the non-detection of millimetre absoiption towards 0132-097 
could be a geometrical effect, where a molecular cloud has a much 
larger chance of occulting the lower frequencies, since the decime- 
tre emission region is likely to be much larger than the millimetre 
region (as may also be exhibited by the wider decimetre profiles, 
Fig-H}- That is, the narrower millimetre emission is less likely to 
intercept the CO and HCO + gas which will be localised in dense 
molecular cores. The fact that the millimetre and decimetre emis- 
sion, and thus absorption, trace different components of the gas 
could be responsible for the significance of the correlation with 
colour not changing when normalising the millimetre column den- 
sities by that of the H I (Sect.l3~2l. 

In the case of 0132-097 , the possibility that the core of the ab- 
sorbing molecular gas is offset from the centre of the emission may 
also result in a reduction of the decimetre covering factor, which 
could be the cause of the normalised OH line strength in 0132-097 
not being as high as expected from its colour (i.e. below the fit de- 
fined by the others in Fig.Q] right), although more absorbers would 
be required to verify this. 

4.2 Location of the absorbing gas 

From the fact that, of our sample, H I has only been detected in 
0500+0 19^ we may be incorrect in our assumption that the host 
galaxy is generally the location of the quasar reddening, as it is for 
only two of the five known OH absorbers. As is the case for the 
remaining three of the known systems, the absorption (and thus the 
cause of the reddening) could be occuring anywhere between us 
and the emission redshift of the quasar. Unfortunately, full range 
spectral scans are not feasible towards most sources with current 
instruments and, over the ranges which are scanned, there is a trade- 
off between the redshift space which can be covered and the depth 
of the search. Thus our motivation for searching for molecular ab- 
sorption w ithin the host galaxy . Spectral scans hav e, however, been 
attempted jMurphv et al.l2003l ; ICurran et al.l2005h , with 0500+019 
being one of the targets common to ours. We summarise the results 
in Table [3] from which we see that < 30% of the redshift space 
towards 0500+019 has been scanned for HCO + and only 12% for 
CO, with some overlap in these ranges. We therefore cannot rule 
out that the cause of the reddening may be at some other redshift 
towards 0500+019 and indeed for the rest of the sample. 

There is also the possibility that, even if the reddening (and 
hence any absorption) is located at the host, it would only be de- 
tected if the optical redshift is known precisely enough to ensure 
that the observed band covers the correct redshift range. The opti- 
cal spectra of these quasars have spectral resolutions of ~ 5 — 8 A, 
which corresponds to a few hundred km s , which is close to our 
observed bandwidths (Fig. |2). However, most of the quasars ob- 
served here have have strong narrow emission lines, which con- 
strain the redshift to better than ~ 10~ 4 and all are quoted to at 



15 As w ell as a possible de tection towards 1107-187, which at z = 
0.48909 iCurranetal]|2011al) . is not covered by the 76.88-77.12 GHz 
(0.4947 < 2 CO < 0.4994) range observed here (Fig. |2}. 0213-026 , 
1504-166 and 1706+006 are undetected in 21-cm ICurranet all 12011 al) . 
with the remainder of th e sample as yet unsearched (see tables 1 & 2 of 
ICurran & Whitindl2010l) . 



least the t hird decimal place, alth ough uncertainties are usually 
not given l lDrinkwater et al .11 1997b . Referring to the ranges cov- 
ered by our observed bands (Table I3.lt , we see that offsets of 
Az « ±0.002 (at low redshifts) to « ±0.01 (at high redshifts) 
are covered, thus making it unlikely that our observations are gen- 
erally tuned to the wrong frequency. 

4.3 Conditions within the host 

Above we suggest the possibility that molecular absorption is not 
detected since the reddening of the quasar light may not necessarily 
be occuring close to the redshift of the background quasar. Over and 
above this, there is the possibility that a large column of absorbing 
molecular gas cannot be located in the host galaxy on the basis that 
the there is a low fraction of cool neutral gas in the hosts of quasars 
and ra dio galaxies with ul tra-violet luminosities of Luv > 10 23 W 
Hz -1 ( ICurran et al 1 120081) . Whether this is due to excitation of the 
gas or these luminosities selecting gas-poor el lipticals (or indeed 
wheth er these two possibilities are inter-related, ICurran &Whitind 
l201Ct) . this appears to be the case for all redshifted (z cm > 0.1) H I 
21-cm absorption searches. 

In order to investigate whether this could affect our sample, 
we use the photometry of th e searched sources (Table |4} to derive 
the UV luminosities (as per lCurran et al.ll2008l) . Showing these in 
Fig. [6] we see the expected increase in luminosity with redshift 
and note that our selection of the reddest PHFS sources (Sect l2,ll > 
means that nearly all of our targets lie below Luv = 10 23 W Hz -1 
(left panel). Since we are unlikely to detect molecular absoiption 
in an environment not conducive to large columns of cool neutral 
gas, in the right panel of Fig. [6] we show the UV luminosity ver- 
sus the colour depth, which we define as the optical-near-infrared 
colour normalised by the depth of the search. As expected from the 
column density — V — K correlation (Fig. @1, t his gives a vertical 
line with this quantity plotted as the abscissqlj and from this defi- 
nition we see that 20 discrete sources (mostly from this work) have 
been searched sufficiently deeply, of which only 12 are located in 
the bottom right quadrant defined by Luv < 10 23 W Hz -1 and a 
colour depth i, 1. 

As expected, in Fig. [6] (right) there is an anti-correlation be- 
tween the UV luminosity and V — K colour, since the former 
is derived from the observed frame optical photometry. As such, 
there is also the possibility that the calculated values of Luv are 
influenced by the dust extinction towards the source. However, the 
fact that 21-cm absorption is never detected above a critical lumi- 
nosity (or below a critical magnitude/dust extinction), as well as 
the cor relation between the 21-c m absorption strength and V — K 
colour dCunan & Whitindl2010l) . indicates that the extinction oc- 
curs within the host in the case of the 21-cm detections, which fol- 
low the expected 50% d etection rate for Luv ~ 10 23 W Hz -1 
dCurran & Whitindl2010l) . Therefore, whether due to high intrinsic 
UV luminosities or a paucity of dust within the host galaxy, it re- 
mains that 21-cm absorption is not detected where Luv ~ 10 23 W 
Hz -1 and where we do not detect 21-cm we do not expect to detect 
absorption by molecular gas. 

Lastly, as well as a covering factor of / ~ 1, as per the 
four known systems, we have assumed an excitation temperature 

16 Adding the three intervening absorbers gives a colour depth of 1.00 ± 
0.02 for both Aco an d 1000 X A HCO +. The intervening absorbers are 
not shown in the plot as, although the ultra-violet luminosities can also be 
estimated, these are remote from the source. 
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Table 3. Summary of the millimetre-band scans towards 050 0+019 with the Swedish ESO-Submillimetre Telescope jMurphv et alJbOoU) in order of 
increasing frequency bands, as tabulated in|Curran et al.| j2005l) . We also show the results for 0727 — 155 since the emission redshift, z e m, of the quasar is 
covered (by CO 1 — ¥ 2, where we show the fractional redshift range based upon z = 1.547 — 1.591). This is therefore included in our analysis (Fig. [3]. 
The column density limits are calculated from the quoted optical depths per a 18 km s _1 (0500+019) and 70 km s _1 (0727-115) channel (Sect. [3~H . 
with T x = 10 K at z = 0. 
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Figure 6. The ultra-violet (A (S 1216 A) luminosity versus the redshift [left] and "colour depth" [right] for the sources listed in Table |4] The colour depth 
is defined by the fit to the CO column density-colour correlation for the four known systems (Fig. [5] left), i.e. [0.77 (V — K) + 12.5]/ log 10 Nrjo an d 
0.77 (V — K) + 12.5/ log 10 (1000 X N-^qq^ ), with the two known associated systems defining the vertical line at a value of 1.0. In the left panel we show 
all of the sources searched, as V and K magnitudes are not required, nor is the re-ca lculated 3cr r.m.s. noise required to be lower than the flux (as is not the 
case for the PKS 0422+004 searches nor HCO+ in Hydra A. lDrinkwater et al.|[l996T) . Note that, after the X 1000 normalisation, the HCO+ points (shown in 
grey for clarity) are located close to the CO points for the two known associated systems. 



of T x « 10 K (at z = 0) [Sect. l4.ll . Since there is no apparent cor- 
relation between the optical-near-infrared colour and the excitation 
temperature for the four known systems (Table IzfH. we cannot es- 
timate temperatures for any of the searched sample. If the gas were 
at higher excitation temperatures, being mostly observations of the 
lower rotational transitions (Table |3Tl . the limits are less sensitive 
than quoted since a larger total column density would be required to 
give the same observed flux in these transitions. This has the effect 
of worsening the limits to the point where only two of the searched 
sources are (just) potentially detectable when T x w 100 K (Fig.[7]l, 
although if the gas is diffuse, such excitation temperatures may not 
be attainable, even in the presence of high kinetic temperatures. 



17 For 1413+135, V-K = 4.81 and T x = 8 K, similar to the excitation 
temperature of 1830-211 with V - K = 6.25. 



5 SUMMARY 



We have undertaken a survey for molecular absorption in the 
millimetre-band at high r e dshift. Li ke all previous surveys (since 
IWiklind & Combeslll995l Il996allbl 1 19971) . we have not detected 
millimetre-band spectral line absorption in any of the targets. In 
the case of absorption due to intervening sources, it has been es- 
tablished that the paucity of detections is due to the traditional op- 
tical selection of these objects, in which the redshift of the inter- 
vening object is usually known, bia sing against the re ddest, and 
thus dustiest, absorption complexes dCurran et al.ll2006t) . We have 
therefore aimed to circumvent this bias, by searching for molecu- 
lar absorption within the hosts of particularly red objects, although 
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Table 4. The magnitudes and calculated A PS 1216 A luminosities [W Hz 1 ] of the quasars and radio galaxies searched in millimetre- 
wave absorption. The bottom panel lists the two known associated absorbers. 



Quasar 




Molecules 


Ref 


B 


V 


K 


Refs 


1o£i n L/JJ\r 
ol(J J - J U V 


4C +40.01 


0.255 


CO, HCO+ 


W96 


17.89 


17.67 


14.87 


Z04.S06 


22.34 


PKS 0113-118 


0.672 


CO, HCO+ 


D96 


18.27 


17.95 


14.37 


Z04 


23.11 


PKS 0213-026 


1.178 


CO, HCO+ 


Cll 


21.33 


20.82 


15.17 


F00 


22.12 


[HB89] 0234+285 


1.213 


CO 


W95 


19.20 


18.35 


13.33 


A09, S06 


23.37 


PKS 0335-122 


3.442 


HCN 


D96 


21.02 


20.11 


17.51 


E05.H01 


23.92 


[HB89] 0422+004 


0.310 


CO, HCO+,HCN 


D96 


13.89 


14.12 


11.31 


F04 


24.55 


PKS 0434-188 


2.702 


HCO+.HCN 


C08 


19.25 


— 


16.24 


E05 


24.14 


PKS 0438-436 


2.852 


CO 


D96 


20.74 


19.91 


16.09 


B04,E05 


23.77 


PKS 0446+11 


1.207 


CO 


W96 


21.43 


20.36 


15.37 


A09,S06 


22.17 


PKS 0454+06 


0.405 


CO 


Cll 


19.38 


18.74 


14.42 


F00 


21.57 


PKS 0500+019 


0.58457 


HCO+ 


Cll 


22.50 


21.35 


15.6 


D97,C03,S96b 


20.38 


PKS 0521-365 


0.0552 


HCO+ 


D96 


15.60 


14.60 


11.33 


H07.S06 


21.96 


PKS 0528+134 


2.065 


CO 


W95 


20.09 


20.00 


15.06 


Z04.G00 


23.51 


[HB89] 0537^141 


0.894 


CO, HCO+,HCN,CS 


D96 


17.93 


17.34 


13.02 


F00 


23.28 


PKS 0601-17 


2.711 


HCO+.HCN 


C08 


20.45 


— 


— 


H01 


23.66 


CGRaBSJ0650+6001 


0.455 


CO 


W96 


20.85 


— 


14.63 


Z04,S06 


21.03 


PKS 0727-115 


1.591 


CO 


M03 


21.13 


20.07 


14.52 


A09,C03 


22.42 


[HB89] 0735+178 


0.424 


CO 


W95 


16.30 


15.68 


13.03 


A09.S06 


23.17 


PKS 0823-223 


0.9103 


HCO+ 


W96 


— 


16.11 


12.34 


F93 


23.51 


Hydra A 


0.0538 


CO, HCO+ 


D96 


13.90 


12.87 


10.90 


V91,C03 


23.89 


[HB89] 0954+658 


0.368 


CO, HCO+ 


W96 


17.18 


16.80 


12.43 


P04,S06 


22.92 


PKS 1026-084 


4.276 


CO 


C08 


21.07 


— 


— 


H01 


24.31 


PKS 1107-187 


0.497 


CO 


Cll 


22.44 


21.10 


15.95 


F00 


19.16 


FTMJ1341+3301 


1.720 


HCO+ 


Cll 


23.39 


21.68 


14.49 


A09.G07 


22.36 


SBS 1347+539B 


0.978 


HCO+ 


W96 


17.95 


17.52 


15.07 


A09,S06 


23.72 


[HB89] 1418+546 


0.152 


CO 


W96 


16.25 


15.60 


11.43 


A09.S06 


22.18 


PKS 1430-155 


1.573 


CO 


Cll 


22.50 


23.24 


17.50 


D97,F00 


21.79 


3C 309.1 


0.905 


CO 


W96 


16.80 


16.57 


14.72 


Z04.S06 


24.27 


PKS 1504-167 


0.876 


HCO+ 


Cll 


20.28 


19.75 


14.01 


F00 


22.36 


PKS 1548+056 


1.422 


CO 


D96 


19.37 


18.65 


14.30 


F00 


23.21 


PKS 1555+001 


1.77 


CO, HCO+ 


D96 


20.34 


19.95 


16.24 


F00 


23.33 


PKS 1622-253 


0.786 


HCO+ 


W96 




20.6 


14.86 


d94,S06 


22.12 


PKS 1622-29 


0.815 


HCO+ 


W96 




18.38 


14.15 


R02,S06 


22.83 


PKS 1725+044 


0.296 


CO, HCO+ 


D96 


17.92 


17.46 


14.13 


F00 


22.06 


[HB89] 1749+096 


0.322 


CO, HCO+ 


W96 


17.98 


17.39 


12.00 


O09,S06 


21.87 


[HB89] 1823+568 


0.664 


CO, HCO+ 


W96 


17.41 


17.02 


13.84 


Z04.S06 


23.18 


Cygnus A 


0.0561 


CO, HCO+ 


D96 


17.04 


15.52 


10.28 


V91.C03 


20.78 


CGRaB J2022+6136 


0.228 


HCO+ 


W96 


19.55 


17.97 


14.28 


Z04,S06 


20.19 


[HB89] 2200+420 


0.688 


CO, HCO+ 


W96 


13.25 


13.32 


10.49 


F04,S06 


25.34 


PKS 2223-052 


1.404 


CO 


D96 


18.59 


18.33 


14.69 


F00 


23.69 


3C 454.3 


0.859 


HCO+ 


W95 


16.03 


16.19 


13.06 


Z04.S06 


24.11 


PKS 2329-162 


1.155 


CO 


D96 


20.72 




16.58 


H01,W83 


22.72 


PKS 2337-334 


1.802 


HCO+ 


Cll 


22.93 


21.89 


16.39 


F00 


22.05 



PKS 1413+135 0.246710 CO, HCO+, HCN, HNC W97 21.37 19.74 14.93 A09.S06 19.41 

B3 1504+377 0.67150 CO, HCO+, HCN, HNC W96 23.24 21.82 16.10 A09,S96a 19.81 



Molec ular search references: D96 - IDrinkwater et al.l (fl 996). W95 - Iwiklind & Combes! jl995h . W96 - Iwiklind & Combed ll996bh 

W97 - IWiklind & Combed fl997l). M03 -iMurphv et aljfeool) C08 -fcurran et al.H2008l) . C ll -this paper. 

Photometry r eferenc es: W83 - IWright et all dl983l) , F93 - iFalomo et alj Jl993l) d94 - Idi Serego-Alighieri et al.l jl994l), S96a - 



^o-Atigr 

IStickel et alj ll996al). S96b -Istickel et all ll996bl) , D97 - IDrinkwater et al.l jl997l), F00 - iFrancis et aU feOOofc GOO - 
feOOOt). E01 - [Ellison et~ai] J200lh. H01 - SuperC OSMO S Sky Survey Hamblv et all 1200 ll). R02 - iRomero et al. 

ICodv & Braunl j2003l). B04~-lBrocksopp et alj teOol), F04 -iFiorucci et alj i2004l) . P04 - IPapadakis et alj <2004l), Z0 4 - Izacharias et al 
l2004h . E05 - Ellison et al] J2005t) . S06 - 2MASS jSkrutskie et alj|2006l) . A09 - SDDS DR7 <Abazaiian et alj|2009h . O09 - lOiha etal 



Ghosh et"ai] 
20021). C03 - 



2009), FPC - P. Francis (priv. comm.). 



like IDrinkwater et al.l d 1996h . who also targetted red sources (and 
type-2 AGN0 we find no absoiption. 

Although as ICurran & WhitindfeOlol) have shown, the bulk cool neutral 
gas (as traced by H 1 21-cm absorption) is not located in the obscuring torus, 
invoked by unified schemes of active galactic nuclei. 



A possible reason for some of the non-detections is that only 
12 of the 40 objects searched are known to be located in the bot- 
tom right quadrant of the ultra-violet luminosity-co/owr depth plot 
(Fig.El right), i.e. where L uv $ 10 M W Hz" 1 and [0.77 (V - 
K) + 12.5]/ log 10 Nco w [0.77 ( V - K) + 12.5]/ log 10 (1000 x 
Nuco+) ~ 1- This suggests that only a dozen of the sources may 
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0-> 
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[0.77(V-K) + 12.5]/log 10 N C0 [10 3 x N HC0+ ] 
Figure 7. The colour depth distribution for various excitation temperatures. 



from the much smaller cross-se ction of the molecula r gas, located 
in small and dense regions (e.g. lLiszt&Lucaslll996h . in compari- 
son to that of the atomic gas. This means that only a narrow sight- 
line will occult the quasar, although the DLA may be apparent 
through the more widespread 21 -cm absorption. 

- As well as covering factor effects, if the excitation temper- 
ature of the putative absorbing gas is higher than the assumed 
T x = 10 K (at z = 0), the column density limits are poorer than 
those calculated here. If large enough (T x > 100 K), none of these 
searches, mostly of the lower rotational transitions, have been sen- 
sitive enough to detect CO or HCO + absorption in these sources, 
although the absorbing medium may be too diffuse to reach such 
temperatures through collisional excitation. 



Drinkwater et al. I d 19961) accounted for their non-detections 
by suggesting that the extinction may be occuring outwith the 
host galaxy or that the X-ray flux from the AGN may be photo- 
dissociating the molecules. These possibilities are similar to the 
po ints made above, al t hough by quantifying these via the findings 
of ICurran et al.l f2 006, 2008j), we find that not all of the targets of 
iDrinkwater et al.l d 19961) are sufficiently faint and reddened to have 
been detected, no matter the location of the obscuring material. 
Therefore, the key to finding new redshifted molecular absorbers 
in the millimetre-band is through the selection of the faintest opti- 
cal objects (giving Luv ~ 10 23 W Hz -1 ) and performing spectral 
scans towards these to colour depth limits of > 1, as defined by the 
column densities and colours of the four known absorbers. In order 
to also circumvent the covering factor effect, pilot searches for OH 
in the decimetre band could be undertaken, with the low rest fre- 
quency (1667 MHz) allowing a full spectral s can in only four sepa - 
rate tunings with the Square Kilometre Array dCurran et al.l2004d) . 



be subject to an ultra-violet flux which permits the presence of cool 
neutral gas dCurran et al.ll2008h . while being "sufficiently red" to 
indicate the p resence of a large c olumn of dust along the sight-line 
to the quasar dCurran et al .1120061) . Given that: 

- For these, the reddening could be occuring anywhere along 
the line-of-sight, and not necessarily in the host galaxy, as is the 
case for three of the five known OH absorbers (four of which make 
up all the known redshifted millimetre-band absorbers). The only 
reason we chose to observe the emission r edshift is that i t give s 
a frequency t o whi ch to tune the receiver 

- lMurohvetaHd2003h ; 
ICurran et alj d2005l) previously attempted millimetre-wave spectral 
scans towards very dim objects, but the redshift range scanned is a 
trade-off with the time spent on each frequency, resulting in poor 
optical depth limits over what were limited redshift ranges in only 
a few sources. 

- This fifth absorber, which was detected in OH but not HCO + , 
may be the result of lower millimetre-wave cov ering factors 
in co mparison to those at decimetre wavelengths dCurran et al.l 
120071) . This may be analogous to the effect seen in the opti- 
cally selected DLAs (Sect. [TJ, where, due to the steeper cos- 
mological evol ution in the heavy e lement abundance in the H2- 
bearing DLAs, ICurran et al. I d2004 suggest that these constitute 
a more ho mogeneous class of ob jects than the general popula- 
tion DLA dProchaska et al.l 120031) . That is, H2 absorption may 
only be observed in a lim ited subset of possible sight-lines and 
IZwaan & Prochaskal |2006) suggest that a distinction could arise 
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